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Continuing our analytic computation of the first-order self-force contribution to the “geodetic” 
spin precession frequency of a small spinning body orbiting a large (non-spinning) body we provide 
the exact expressions of the tenth and tenth-and-a-half post-Newtonian terms. We also introduce a 
new approach to the analytic computation of self-force regularization parameters based on a WKB 
analysis of the radial and angular equations satisfied by the metric perturbations. 


I. INTRODUCTION 


The impending prospect of detecting gravitational- 
wave signals from coalescing compact binary systems mo¬ 
tivates renewed studies of the general relativistic dynam¬ 
ics of binary systems made of spinning bodies. It has 
been emphasized in Ref. [lj that a simple way of com¬ 
puting (to linear order in each spin) the spin-dependent 
interaction terms Hi nt = Clf° ■ Si -ffif 0 ■ S 2 in the Hamil¬ 
tonian of a binary system was to compute (when consid¬ 
ering, say, the term linear in the spin precession an¬ 
gular velocity of Si in the gravitational field generated by 
the two masses mi, m 2 , and, eventually, the spin S 2 . In¬ 
deed, this spin precession angular velocity (which can be 
obtained by writing that Si is parallely propagated along 
the world line of mi) is simply equal to the coefficient 
of Si in Hint• On the other hand, it was recently re¬ 
marked that, in the simple case of a binary moving 
on circular orbits, the (^-component of the) spin preces¬ 
sion, fli°, could be expressed in terms of the norm |Vfc| 
of the covariant derivative of the helical Killing vector 
k = dt + Cldc/, characteristic of circular motions, namely 

cif° = n - |vfc|, (1.1) 


where Cl denotes the orbital frequency. [The gauge- 
invariant quantity |VA;| can be viewed as a first- 
derivative-level generalization of Detweiler’s redshift in¬ 
variant [3], which is expressible in terms of the norm |fc| 
of the Killing vector k] 

The gauge-invariant functional relation between Clf°, 
or equivalently |VA;|, and the orbital frequency Cl has 
been recentlvstudied (both numerically and analytically) 
in Refs, jg, Q. In particular, we have derived (as part 
of a sequence of analytical gravitational self-force stud¬ 
ies) in 3 the first-order self-force contribution (linear in 
the mass ratio q = mi/m 2 <C 1) to the “geodetic” spin 
precession frequency Clf° to the eight-and-a-half post- 
Newtonian (PN) order, i.e. up to terms of order y 8 ' 5 
included, where 


V = 


( Gm2Cl 

\ 


2/3 


( 1 . 2 ) 


is a convenient dimensionless frequency parameter of or¬ 
der 0(l/c 2 ). [We henceforth use, for simplicity, units 


where G = c = 1.] As in Q we restrict ourselves here to 
the case of a small spinning body mi, Si, orbiting a large 
non-spinning body m 2 , S 2 = 0. 

The aim of the present note is to report on an ex¬ 
tension of our previous analytical computation of spin 
precession to the 10.5PN level, i.e. up to terms of order 
yio.5 i nc i uc i ec b This extension was motivated by private 
communications from Dolan et al. 5| who pointed out 
apparent discrepancies (starting at level 0(y 7 )) between 
some of their high-accuracy numerical results (see Table 
III in Ref. Q) and our published 8.5PN analytical re¬ 
sults. These discrepancies led us to carefully re-examine 
our previous computations, and to push them to higher 
PN orders. We so discovered that, though all our ba¬ 
sic analytical building blocks were correct, their manip¬ 
ulation by an algebraic software led to some instabili¬ 
ties (due to the length of the analytical expressions at 
high PN orders), that had led to a few errors in our fi¬ 
nal results. More precisely, the rational term, among the 
seven (transcendental) contributions to the coefficient of 
y 7 , was incorrectly obtained, and, in the coefficient of 
y 8 (which contains fifteen different contributions), both 
the rational term and the coefficient of 7r 2 were incor¬ 
rectly obtained. Correspondingly, there were errors in 
the (rational) coefficients of y 7 and y 8 in the subtraction 
term B{y). [See detailed results below.] After having 
found these errors, corrected them, and communicated 
the corrections to Dolan et ah, the latter authors con¬ 
firmed that our 0(y s ) corrected results were now in satis¬ 
factory agreement with their high-accuracy numerical re¬ 
sults. [More recently, Shah 7| independently pointed out 
to us the three discrepant coefficients mentioned above, 
which we had already analytically derived, and which he 
and his collaborators had independently derived by using 
the numerical-analytical method of Ref. Q.] 


II. TECHNICAL REMINDERS 

Let us recall the notation and main technical results of 
Ref. H that we shall need to express our new results. We 
consider a two-body system of masses m\ and m 2 , mov¬ 
ing along circular orbits, in the limit mi C m 2 ■ Here we 
only endow the small mass mi with spin S \, keeping the 
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large mass m 2 non-spinning. This means that one is deal¬ 
ing with linear perturbations h^(x x ) of a Schwarzschild 
background of mass m 2 by a small mass mi, moving on 
a circular orbit of radius tq. As emphasized by Detweiler 
3, the perturbed metric admits the helical Killing vector 
k = dt + kldc/,, i.e., the metric perturbation depends only 
on <j) = <j) — fit, r and 9 , r, 9). 

The four-velocity of mi, normalized with respect to 
the metric x x ) = gffi + qh + 0(q 2 ), (here q = 
m\/mi <C 1 and the superscript R indicates the regular 
part [9] of hfj V (x x ) around the world line of mi), can be 
written as 


bk 

v ? = Wi sTk "’ 


r — — 

~ |*|’ 


( 2 . 1 ) 


where (to linear order in q) 

|fc| = ^[~9%ki i k v ] 1 = ^Jl- - fl 2 rg - qhkk 


= i/l- 


^ - Wr 2 
TO 


1 -\q 


kkk 


2 y 1 _ 2^2 _ ^ 2 r 2 


( 2 . 2 ) 


with h kk = {x)k^k v ]\. Writing that mi moves 

along an equatorial circular geodesic yields the conditions 
d^g kk = 0) which lead to [|| 


n = 



1 ~q 


4m; 


■[drA 


kk Jl I 1 


[d^hf k ] 1 = 0 . 


(2.3) 

(2.4) 


Eq. (12.31) allows one to trade the gauge-dependent ra¬ 
dius ro for the gauge-invariant dimensionless frequency 
parameter y, Eq. & using 


m 2 m 

r 0 = - q 

V 


6 y 3 


[drh kk \ 1 , 

^ = V + q^ldrhf^. (2.5) 

The geodetic spin-orbit precession frequency along the 
world line of mi has, as only nonvanishing component, 
fif° = nf° given by Eq. CUD above. In this equation, 
the norm |Vfc| of the covariant derivative of the helical 
Killing vector k = d t + fid# is defined as 


|v*| 2 = ^(v M M(v' 1 fc ! '), 


( 2 . 6 ) 


where all tensorial operations are done with the metric 
ffuv( x )- The explicit expression of |Vfc| can be written as 


where 


is the well-known result for gyroscopic precession (with 
respect to a rotating, polar-coordinate frame) in a 
Schwarzschild background m , and where 


-,3/2 


y 


htcf) 


y 


m 2 (l - 2 y) 

' ~ Orh^k) • 

y 


2(1 — 2?/)(l — 3y) 


kk 


(2.9) 


In Eq. (12.91) all quantities are to be regularized and eval¬ 
uated for 9 = 7t/2. 

The quantity 5(y), which measures the fractional first 
order self-force (1SF) correction to |V*|, is equivalent to 
the quantity 5i/)(y) which measures the 1SF contribution 
to the dimensionless ratio |2j 

V>(y) = ^ 3y[l +q6(y)+0(q 2 )}. 

Explicitly, we have 


( 2 . 10 ) 


Sr/j(y) = -yj 1 - 3 yS(y). 


( 2 . 11 ) 


Following the methodology explained in Refs. PH3, 
and extending the results of Ref. [3j to higher post- 
Newtonian orders (by using radiative solutions, 

A( up ), up to l = 7), we have computed S(y) up to or¬ 
der y 10 - 5 . 


III. NEW HIGHER POST-NEWTONIAN 
TERMS IN 6{y) AND 5rp(y) 

Before listing the complete expressions of 5(y) and 
5i/j(y) to order y 10 - 5 let us indicate that our previous 
O (y 8 ' 5 )—accurate results missed one term at level y 7 and 
two terms at level y 8 , while the y 7 ' 5 and y 8 ' 5 terms were 
complete. 

More precisely, the correct 0(y 8 5 )—accurate expres¬ 
sion of <5(y) is obtained by adding Acf y 7 + Acf y 8 to Eq. 
(4.33) in ||, where 


Acy = — 

A Co = 


1485630311863 

45831035250 ’ 

8 0 25377697082469367 


8505 


262980313505910 


(3.1) 


Equivalently, the correct 0(y 8 5 )—accurate expression of 
6ip(y) is obtained by adding Ac^y 7 + Ac^y 8 to Eq. (5.4) 
in 0, where A c 7 ^ = —A c? and A = — Ac| + | A cf, 
i.e., 


|Vfc| = |Vfc|(°)(l + y<5(y) + 0(y 2 )), 

(2.7) 

Ac) 

|V*|® =9y/l^%, 

(2.8) 

Ac! 


_ 1485630311863 


45831035250 




8505 


7G + 


629539392522290711 

13149015675295500 


■ (3-2) 
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The full 0(y 10 ' 5 )—accurate expressions of S(y) and 5ip(y) read 


$(y) = 


0 


15462423 172 

-1-7 - 

4480 5 1 


X 




30832 

^oT 7 ' 


26536 10/0 

- ir i/z tt 

1575 y 

/16156122817 
V 1209600 

670667 ir/o 

- 

22050 

/ 41432062371919 96697099_ 58208 2 1291394011 1007542476707 2 

^ 2540160000 + 141750 7_ 105 ln( ' 3638250 353894400 ^ 

9765625, . , 96697099, . . 2364633, 856, . , 2 162286431837 4 869696, 

^^ ln(5) + ^8350^ ln(2/) + ^232^ ln(3) “^ ln(y) + 335544320 * ~^57T ln(2)7 

iQ/i/i q 7,0/1 x 3424 2 434848, , .A 8 


V 2540160000 141750 ' 105 3i 

9765625, . , 96697099, . . 2364633, 856 

~ 28512 ~ ln(5) + ^8350T ln(y) + 1232^ ln(3) " IT 
1344 3424 , , . 3424 2 434848, . ,, . / 

^C(3) - - 25-7ln(y) - ^7 2 - ^ M2) W 

0 070 rr A 0070 


3344, . , 512537515 2 15416, , , 40581, 1407987 A 

-In 2)- tt 2 -In (y) -In 3) +- tt 4 - 

21 v ; 393216 105 140 v ; 524288 J ‘ 

1007542476707_ 2 
353894400 ^ 

i no7 n n 


O AO 

3872542979 17/2 

- y 1 7r 

13097700 ' 


2/° 


+ 



7776 

+ 

5 < 


32288 


75 


28431 


49 


118580138377_ 45728 2 58794404629417 

’ 14553000 7 + 1225 n( ^ + 3972969000 ^ ) 

71 2 20486328125, , r , 118580138377, . , 

_7r - 5189184 ln(5)+ 29106000 Hy) 



(3-3) 
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87055 , 3772 
28 + T05- 


y KJKJ 

3772, . , 

w ln < 2 >- 


26536 ,o/9 

- V ld ' Z TT 

1575 y 

( 149628163 7628 3814, 

h { -S900“ + — 7+ —' 

113411 / 

22050 y 

( 74909462 58208, 

V 70875 ' 105 w 

160934764317 4 1344 8 

335544320 ^ 5~'^ + “ 

9765625, /r , 856, . , 2 43^ 

-ln(5) H- My) H- 

28512 w 25 y ’ 15 

403109158099 3424 , . A , 

9922500 25 ' yy> J y 

1179591206 ,,,, 


7628 3814, . , 4556, 12879, 1407987 4 297761947 2 \ 

IT 7 + If ln(2/) + IT hl(2) + "ST" ln(3) - ”524288" 71-4 + ^9321^"') 3 

58208, , 9 340681718, 164673979457 2 

-ln (2 2 H-In 2 ) +- 7 r 2 

105 v ' 1819125 v ; 353894400 

1344 869696, 3424 2 199989, 

-C(3) H-In 2 7 +- 7 2 -n 3) 

5 1575 w ' 25 1 352 y ’ 

556, . 434848, . . 37454731, , . 

+ 

3424 , , A . 


yyz/iouu zo j 

1179591206 17/2 

H- y 1 ' 2 tt 

3274425 y 

f 4454779894_ 1064368 2 138895624334 22832200546571 2 

+ V 606375 7 1225 Uy ’ 9029475 n( ) 8808038400 n 

758053590944149 4 3077728, 3376 2 28431, /ox2 

- 7 r 4 - 1152C 3)-ln (2 7 H- 7 2 H-ln(3 2 

21474836480 sv ; 11025 w ’ 15 1 49 w 

71602663581, 11576171875, , r . 844, . x2 28431, . 

-ln(3) H-ln(5) 4 -ln(w ) 2 H-ln(j/ ln(3) 

7847840 v ' 2594592 y J 15 yy 49 yy y 1 

1538864, . . . 2227389947, , , 56862 , . , 3985926908910281 

-ln(2) ln(w)-ln(w) H- 7 ln( 3 ) H- 

11025 v yy ’ 606375 yy> 49 ' w 1146048750 

56862, , . 3376 , . A g 

+ —Jg- ln (2) ln (3) + -^-7 In {y)J y 

( 660044682996077 434848 , 46528736 6228256, 

v 524431908000 4725 165375 ' 11025 v ’ 

23264368 , . A 19/2 

/11467229058074 21138410295 30719079112, 1306135539288758, 

+ - 7 H- 7 T 6 +-In 2) 2 -In 2 ) 

V 496621125 134217728 3274425 ' ; 147496474125 v ; 

152033994681460553 2 755954175166870909 4 680336,478423984, . . 

- 7 r 2 -7T 4 +-C(3-In 2 7 

13317754060800 2748779069440 105 sw 654885 yn 

35570296 2 54832464, . , 2 214411899501351, 437134765625, /r , 

11025 7 - ln(3) + 3453049(500 ln(3) “ 2075(5736 ln(5) 

—TToST■"« -1SmT ln ^ hl<3 > - -65i«r hl(2 > + 496621125 

109664928 , 552424223705497767347 109664928, 35570296 , . , 

-13475 -7 ‘" (3) + 20649506377500 -1347^ ‘" <2) ‘" <3) “ TToffi -7 

678223072849, .A in 
~ 370656000 ‘" (7) ) ^ 


” 370656000 M(7) K 

/ 178279193702345741_ 11255086 

_ ^ 26512946460000 33075 

369603 , 23162039 , . A 21 

-7rln(3) H-7rln(j/) ] y zl 

343 ' ; 128625 yy, J y 


46324078 2889622 , . . 

-7T7-In 2)7r 

128625 ’ 7875 v ; 

oUv 11 ) ■ 
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Note that, in the above expressions, we used the 
computer-algebra-related notation ln(a)" to denote 
In" (a) and 0\ n (y n ) to denote a term of order y 11 In" y for 
some n. The corresponding 0\ n (rt 11 )—accurate expan¬ 
sion of the effective gyrogravitomagnetic ration gg SF (u) 
is given in the Appendix. 


IV. ANALYTIC EXPRESSION OF THE 
SUBTRACTION TERM 

Prompted by Dolan et al. 5], who pointed out dis¬ 
crepancies at order y 7 and y 8 between our Eq. (4.30) 
in |3| and their (unpublished) corresponding expression 
for the subtraction term B(y), we have found a way to 
derive an exact analytic expression for B(y) within our 
formalism, which is based on Regge-Wheeler-Zerilli-type 
tensorial multipolar expansions. As we shall now ex¬ 
plain, our derivation is a novel approach grounded on 
a Wentzel-Kramers-Brillouin (WKB) analysis of the ho¬ 
mogeneous radial (Regge-Wheeler, RW) equation satis¬ 
fied by the fundamental building blocks, A'; n and A up , of 
our formalism. This WKB approach (which we explain 
in detail below) is quite different from the approach tra¬ 
ditionally used in gravitational self-force theory, which is 
based on local, Hadamard-type expansions of the metric 
huv, in Lorenz-gauge, near the world line of m\ (see e.g., 
[3, [la, [l3| )■ In addition, our approach defines the sub¬ 
traction terms by considering the limit l —> 00 where l 
denotes the degree in a tensorial multipolar expansion, 
while the usual self-force calculations define subtraction 
terms by considering a limit l s —> 00 , where l s denotes 
the order in a scalar multipolar expansion. One can show 
that, for the quantities we shall consider, the two differ¬ 
ent limiting procedures should give the same subtraction 
terms at leading order. [However, at higher orders in lo¬ 
cal singularity expansions, the extension ambiguities of 
such expansions do not imply anymore their equivalence.] 
Let us start by recalling the form of the WKB approxi¬ 
mation of the solutions of a one-dimensional Schrodinger 
equation, say 

= (4.1) 

The WKB solutions of Eq. (14.11) are written in the form 

'R(rc) = e %+ s i+ 0 ( s ). (4.2) 

As indicated here, it will be sufficient for our purpose to 
keep only the leading and next-to-leading terms in the 
WKB expansion. At this order of approximation, the 
two independent solutions of Eq. m read 

*±(x) = C ± e±I / ^ X , p(i) = V / QW, (4-3) 

Vp( x ) 

corresponding to 

Y =± J P( x ) dx - Si = -iln p(x). (4.4) 


The choice C± = 1/ \/2 would imply that the Wronskian 
of these solutions is 1 : 


W = ’L_’L+ - = 1. (4.5) 


Note that we will use the WKB approximation in the 
classically forbidden domain, where Q(x) is positive so 
that the solutions T-t are exponentially growing or de¬ 
caying. 

We first apply this approximation to the (homoge¬ 
neous) radial RW equation 


dr * 2 


1(1 + 1 ) 


x = m{^-^)-v 2 rn*ny. 

( 4 - 6 ) 

Here, 77 = 1/c, m is the spherical harmonics order, Q 
denotes the orbital frequency, and 


dr* 


dr 

/O’) ’ 


/O’) = 1 - 


2 Mr/ 2 
r 


(4.7) 


The spatial variable (denoted x in Eq. (14.11) 1 in this one¬ 
dimensional Schrodinger equation is r*, while we shall 
take as small expansion parameter h the quantity 


h 


\ 

L’ 


(4.8) 


where we introduced the convenient notation 

L = l + \ • (4.9) 

Note indeed that the coefficient 1(1 + 1) in the centrifugal 
potential can be written as 

l(l + l) = L 2 - i, (4.10) 


and is of order ~ p. 

In order to capture the near-world-line singularity ex¬ 
pansion within our tensorial multipolar expansion, we 
need to consider a limit where both l ~ L and m tend 
to infinity with the ratio w = m/L being kept fixed. In 
this limit the two dominant terms (of order 1 /h 2 ) in Eq. 
(14.11) are 


7? = ,(, + 1) 


3 -/( r ) - V 2 


1(1 + 1 ) 


fr 


O(L 0 


= L 


= L 2 


m 




1 


f(r) - r] 2 w 2 n 2 


O(L 0 ) 


+ O(L 0 ). 


(4.11) 


Correspondingly to the accuracy used in Eq. (14.21) . we 
can neglect the terms of order O(L 0 ) in the above equa¬ 
tion, which notably means neglecting the term 6 M /r 3 in 
Eq. (14.61) . At this stage no expansion is performed in the 
PN-parameter rj = 1/c. 

Introducing the notation 


A(r) = 1 — 


2Mif 

r 


2 2n2 2 

rfw irr . 


(4.12) 
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we have 


Q* = vl = L 


2 A 0) 


p* = L 


V A ( r ) 


(4.13) 


so that the building blocks of the WKB solution (14.31) 
read 


y/A(r) dr 
f(r) r 

(4.14) 


(4.15) 


(4.16) 


(4.17) 

where we have re-absorbed the factor yfh in the constant 

C± [c ± = c±/Vl\. 

The final result of this WKB analysis is that two inde¬ 
pendent solutions of the RW equation (14.61) are 

X±=C ± j£ I e ±L f^ . (4.18) 

We checked that the PN expanded solutions of the RW 
equation that we constructed in our formalism (Ml 
agree with those WKB solutions, with the following cor¬ 
respondence 


A V) R 

„ vZ r L fY*dr 

- A 1/4 C 

A'( U p) ?= 

s -Lj^dr 

AV4 e 


Note for instance that, when expanding in powers of y 
the right-hand-side (rhs) of X( in ), as given in Eq. (14.191) , 
its leading order is y/re Llnr = r l+1 in agreement with 
the normalization of our PN solution which was chosen 
as 


= ± [ p*dr* = ± [P*~r = ±L 


and 


More explicitly 


so that 


Si = -- lnp* . 


Vp: = VZ 


C± 


= C± 


A 1 / 4 


y/K M/4 


4“n)( r ) =ri+1 ( 1 + ^“( r )) > ( 4 - 20 ) 

with A lu (r) = 0(y 2 ). 

Inserting the above WKB solutions for X( in ) and W up ) 

in the analytical expressions for ^±( odd / evcn ) gj ven j n R e f. 
3 [see Eqs. (4.10) and (4.12) together with Eqs. (4.11), 
(4.13) and (4.23) there] yields expressions for 5±(y) = 
S±’ even (y) + 5±° dd (y) of the form of Eq. (4.28) there, i.e., 

s Hv) = ^2 s L(y) = ±LA (y) + B( ,y ) + ■ 

m ^ ' 

(4.21) 


At this stage the subtraction term B{y) is given by a sum 
over to of the form 


E/(?)i«.i’+»(?) 


d^ lm 

de 


(4.22) 


where Yi m (9, <j>) and its 6— derivative are both evalu¬ 
ated at 9 = 7t/ 2 (and (j> = 0). Such a sum can be 
asymptotically evaluated, in the limit L —> oo with m/L 
fixed, in terms of an integral, between —1 and 1, over 
the variable w = m/L. In order to do so one needs 
asymptotic estimates for |Yz TO | 2 and as functions 

of w in the large L limit. Such asymptotic estimates 
can be derived by a WKB analysis of the 9 differen¬ 
tial equation satisfied by <di m (9) (defined by factoring 
Yi m (9,(j>) = 0z m{9)e im ^). Indeed, 0/ m (0) satisfies a one- 
dimensional Schrodinger equation of the type ED, when 
using the variable A = J , 2 d9/ sva.9, namely 

-^&lm = ~P 2 WQlm, (4.23) 

with 


P 2 (A) = 1(1 + 1) sin 2 9(A) — m 2 . 
This leads to WKB solutions of the type 


©Zm(A) = C+- 


z ifP(\)dX 


C-- 


0 -if P(\)d\ 

VpW 


(4.24) 


(4.25) 


for appropriate choices of the constants C± determined 
by regularity conditions at A = — oo (corresponding to 
9 = 0) and A = +oo (corresponding to 9 = tt). When 
evaluating 0/ m (A) and ^0/ m (A) = sin 9-^Q lm at 9 = 
7t/ 2 (i.e., A = 0) one finds the following WKB estimates 


\Yi m & ,0)| 2 « ±-=L=S?%, (4.26) 

2 y/l - w z 

and 

\deY lm £, 0)| 2 « ^y/T^Pd^t- (4-27) 

Z 7T Z 

Here 6^2 (^°-„) is equal to 1 when l — m is even (odd) 
and to 0 otherwise. The above estimates are not a priori 
uniformly valid in the full range —l<m<l because our 
WKB analysis requires L 2 — to 2 1. However, they can 
correctly evaluate the asymptotic values of the integrals 
that we shall be interested in below (which have only a 
relatively small contribution from the neighborhoods of 
the boundary points iu = ±1). We have checked the esti¬ 
mates (14.261) and (14.271) by using the explicit expressions 
of Km(§ .0) and dgYi m ( |,0) gi ven in Eqs. (32) and (33) 
of [H] . A consequence of Eqs. (14.261) and (14.271) is that 




v- 4 tt \d e Y lm (%,0)\ 2 

-7,- 9 


21 + 1 


L 2 


(?) 


1 
7r 


dw y/l — w 2 g (w) 


(4.28) 
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As an example of the application of these asymptotic es¬ 
timates we have computed the analytic expression of the 
L —> oo limit of the first-order self-force redshift quantity 
hkk- Starting from Eqs. (29) and (30) of Ref. |lT| one 
finds that B hkk = lim^ 00 (/iS + h kk,tl) is S iven b Y 


tegral of the first kind (with w = sin a): 


EllipticK(fc) = 


7r/2 


da 


Vl - k 2 


sin 2 a 


(4.31) 


Bh kk 


2y(l — 3j/) 3 / 2 1 I' 1 dw 

\/l — 2 y 7r J _ 1 y^(l — w 2 )(l — k 2 w 2 ) 

22/(1 - 3 j /) 3 / 2 2 


V 1 - 2y tt 


EllipticK(fc). 


where 


fc 2 = 


2/ 


1-22/ 


(4.29) 


(4.30) 


and where EllipticK(fc) denotes the complete elliptic in- 


This result agrees with the subtraction term obtained by 
the usual self-force Hadamard-type analysis I, 03, i-e-, 
the term denoted Dq = (l — 3y)D Q in ED, EH- [Note that 
there is a misprint in the last term of Eq. (56) in [T3 : 
the coefficient of u 7 should read +4409649/524288]. 

When applying the above WKB asymptotic estimates 
(for both the radial functions X±(r) and the angular 
functions Y) m and dgYi m ) to the l —> oo limit of the 
quantity S±(y), (14.211) . we obtain the following analytic 
expression for the O(L 0 ) subtraction term B(y) 


%kb {y) = -J ] t ? [(4 - 92/)EllipticE(fc) - 2(2 - 52/) EllipticK (fc) ] . 

7T V 1 — 22/ 


(4.32) 


Here EllipticE(fc) denotes the complete elliptic integral 
of the second kind 


/.tt/2 - 

EllipticE(fc) = da \/1 — k 2 sin 2 a . (4.33) 

J o 


The expansion in powers of y of Bwkb reads, up to 
the 11 PN level 


Bwkb(u) 


2 4" 128" 1024' 

126849 16567767 7 

- V b 4 - y 7 + 

32768 9flQ71K9 y 


oa/uo 2097152 V 335 
77104836855 9 , 35027350019' 

2147483648 V + 4294967296 

26812467118879 n 

- y 11 + 0(y 

137438953472 y 


,, 12 \ 


In our previous work Q the subtraction term B was not 
derived independently of our computation of but was 
obtained from the large l limit of the PN expanded ver¬ 
sion of ■ The algebraic-manipulation errors mentioned 
above induced corresponding errors in our previous eval¬ 
uation of the PN expansion of B (as pointed out to us by 
Dolan et al. !)• More explicitly, the coefficients of y 7 
and 2/ 8 in Eq. (4.30) in [1} were in error, and Eq. (14.341) 
gives instead their correct values. 


V. CONCLUDING REMARKS 

The analytic computation of the post-Newtonian ex¬ 
pansion of the first-order self-force contribution to spin 
precession has been raised here to the ten and ten-and- 
half post-Newtonian level. 


Our analysis has also corrected two terms (at the PN 
levels 7 and 8) among our previous 8.5 PN-accurate 
calculation of spin-orbit effects [3j. More precisely, we 
have shown that Eq. (4.33) in Q needs to be aug¬ 
mented by the two terms in Eq. (ED. Equivalently, 
Eq. (5.4) in [3] needs to be augmented by the two terms 
in Eq. m • These missing terms were caused by al¬ 
gebraic errors in the manipulation of large analytic ex¬ 
pressions. Note that these errors affected only a few 
terms among many contributions (essentially only ra¬ 
tional terms). The missing contributions to the coeffi¬ 
cients of y 7 and y 8 in Stf>(y) are numerically equal to 
Ac 7 = 32.41537757 and A cf = 47.86801827. These val¬ 
ues are rather small compared to the corresponding typ¬ 
ical values of the general PN coefficient c)([ ~ —0.12 x 3 n 
(linked to the pole singularity of 5 ip{y) at y = \, see Eq. 
(5.14) in 0). The fractional modifications brought to 
the coefficients gl and in Eqs. (6.36), (6.37) and (6.38) 
are correspondingly small, Sge/ge ~ —0.02487821950 and 
Sgr/g 7 ~ 0.0001739775786. As a consequence, correcting 
these terms does not affect any of the significant con¬ 
clusions we reached in ;3| which were mainly aimed at 
describing strong field effects. In particular, our fits Eqs. 
(5.11) and (6.39) did not make any use of the y 7 and y 8 
coefficients but only relied on 3PN information and on 
the strong held numerical data of [2]. 

Finally, we have introduced here a new method for an¬ 
alytically computing the subtraction terms of self-force 
quantities. Instead of the traditionally used Hadamard- 
like near-world-line singularity expansions, our new 
method is based on a WKB analysis of both the radial 
and angular equations satisfied by the metric, when con¬ 
sidering them in the limit l —> oo with the ratio m/l fixed. 
We have shown on two examples (hkk and the spin pre- 
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cession) that our method leads rather simply to closed Appendix A: Higher PN terms in gf: 

form expressions for the subtraction terms. 

Combining the O(j/ 10 - 5 )-accurate computation of S(y) 
above with our recent O(y 10 5 )-accurate computation of 
the main effective one-body radial potential a(u ) fl5j . we 
can raise the PN expansion order of the effective gyro- 
magnetic ratio gg SF from the 0(u 7 5 ) level given in Eq. 
(6.37) of Q to the 0(u 9 ' 5 ) level. We list below the final 
result, expressed in the effective one-body radial variable 
u. 


9 


9s7(u) = 


3 39 2 

-~u - —u 

4 4 


41 

32' 


^ 2 - 


7627\ 


192 


+ —487 + 


23663 


2048 


7T 2 — 


1456 1017 \ 4 

ln ( 2 ) “ ~^FT ~ 24 M«) u 


15 


-fM3) 


9832 


35 


-7 + 


712905 
4096 7 


20 
70696 
+ 105 


M2) - 


161160813 4916 


89600 


35 


ln(u) 


- 7 TU 


93304 
1575 
/315657 

V 280 

4596019 


11/2 
ln(3) H 

13/2 


480829_ 
2835 


16790137 4 674904611 2 2954531 

-7T- Tt" — 


1048576 


7077888 


2835 


M2) - 


18167439833 480829 


7257600 


+ 


5670 


1h(m) 


+ - 


- 7 TU 

12600 

12227517 


3080 
1167584 
525 
451802734 
121275 
118299749 


ln(3) - 1088((3) - 


1953125 

3564 


ln(5) — 


903605468 
121275 


-7 + 


58208^ 

105 


ln(2) 2 - 
ln(it)- 

15/2 


7532631301 


9175040 
58208 


499904, , , 5587843424, , , 

-M2)7- ln ( 2 ) + 


225 


779625 


204902966117_ 4 
335544320 ^ 
48146264595158227 
625895424000 


105 


, , . 249952, , . 14552, , , 2 

7 ln(u) 4 -——— ln( 2 ) ln(u) H— ln(w ) 2 


225 


105 


- 7 TU 

2182950 

52964727700527 

3139136000 

1135089788764019 

42949672960 
164036944, . . 

M2)7 ■ 


hi(3) 


7T 4 — 


141648 

20022888 


C(3) 


366384765625 

41513472 


M5) 


-7 


11025 
4189028005087 


, . , 2 142155 , 

1225 hl(2) + ^“ 7ln(3) + 

90305230479881 ln(2)+ 142155 ln(3)2 


4204284206047 
264864600 

142155 ln(2)ln(3) 


10974904 , 

3675 ' 

1241427590810221 


49 w w 369937612800 
1094977266529990589159 


3972969000 


98 


529729200 
/142517152 


, . , 10974904 , . , 82018472, , . , . , 

ln(u)-——— 7 ln(u) — ln( 2 ) ln(u) 


427173626880000 
142155 


3675 


11025 


98 


2743726, , , 2 
ln(u) ln(3)-——— ln(it ) 2 


3675 


+ 


55125 
/15312301495292259 


, 6965217870900563 213592544 

In(2)7r--— - 7 r 4- 777 ■ 


69060992000 


M3) 


762810048000 
1294640 


165375 


1996192 

4725 


106796272 

165375 


7rln(u) ) u 17 ' 2 


81 


13759767578125, , r , 16569352454284793 

C(3)- ^ nnonooo ln(5) + onooono ^ oonn 7 


4477353976 
1403325 


691974898583334793 


249080832 
540743945464 


-7 


5497558138880 


1368675122796890401 


"77 + 


264895105264 


146495294668800 9823275 

78318056677502249, 847418031, /0 , 2 

■ ln(2)-ln(3 2 

V 1 53900 V ' 


ln( 2)7 - 


9823275 

207604582525 


202280878800 

, , , 2 847418031 , 847418031, 

ln(2 ) 2 - 7 ln(3)- ln(2) ln(3) 


26950 


26950 


3003559322617, 

7 T_ _ In (7) 


402653184 1111968000 

849725095980588589949507303 16550580788732153 


66987659472814080000 


7079830758000 

4477353976 , , , 132447552632, . . , , , 847418031, 1119338494, , x2 

7 ln(u) H-ln(2) ln(u)-ln(u) ln(3)-ln(u ) 2 


404561757600 
9 


ln(w) 


1403325 
78150479 


4410 
1137638861 
308700 


ln(2)7 


9823275 
298035034972802327 
11783531760000 


53900 

1137638861 

7T- T~Z~TTTZ~T -7T7 ' 


1403325 


u 


154350 


4065633 
1372 7 


ln(3)+ 


23679911 
13230 7 


7 r ln(u)J u 19/2 + O ln (u 10 ). 


(Al) 


Here the coefficients of u e and u 7 differ from the ones 
given in Ref. [3| because of the corresponding change in 
Sip. More precisely, the terms to be added to Eqs. (6.36)i 
and (6.37)i in 0 read 


Acf* = Ac d r 


,541 


(A2) 


and 


Acf* = Acf" - |Acf 


(A3) 
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Let us also point out a misprint in the expression of g 7 11 
given in Eq. (6.37 ) 2 of Q: the additional term 


+ 


249952 

225 


In 2 , 


(A4) 


was inadvertently omitted. 
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